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1.  INTRODUCTION 


Recently  developed  numerical  methods  have  been  applied  successfully  to 
performance  prediction  of  antennas  in  complex  scattering  environments.  The 
general  framework  for  these  techniques,  at  frequencies  for  which  the 
antenna*  is  at  most  a few  wavelengths , is  the  method  of  moments  (moment 
method)  as  set  forth  by  Harrington  [Harrington  1968] . Many  computer  programs 
have  been  developed,  based  on  the  moment  method.  A number  of  them  are  appli- 
cable to  solving  antenna  problems  of  quite  arbitrary  geometries.  In  particu- 
lar the  programs  based  upon  thin-wire  modeling  are  applicable.  Not  all 
antennas  are  collections  of  thin  wires,  however.  Hence,  there  are  two  con- 
cerns regarding  these  methods:  (a)  the  accuracy  of  the  model,  and  (b)  the 
accuracy  of  the  analysis  of  the  model.  An  accurate  model  is  one  which  ex- 
hibits all  important  aspects  of  the  actual  antenna.  This  accuracy  can  be 
gauged  either  by  (1)  comparison  with  experiment  or  other  analysis  if  results 
from  these  exist,  or  (2)  continuously  increasing  the  complexity  of  the  model 
(add  wires,  for  example)  until  computed  parameters  such  as  driving  point 
impedance,  radiation  pattern,  or  current  show  little  change.  Gauge  (2)  is 
a convergence  test. 

An  accurate  model  analysis,  on  the  other  hand,  is  not  related  to  the 
relevance  of  the  model,  but  rather  to  the  accuracy  with  which  an  analytical 
tool  predicts  important  model  parameters.  For  example,  a thin-wire  moment 


The  word  "antenna"  as  used  throughout  this  report  is  intended  to  mean 
"antenna  and  its  environment". 
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method  computer  program  could  be  applied  to  a wire-gridded  model  of  a 
corner  reflector  ar.tenna  to  determine  its  radiation  pattern.  The  result 
might  accurately  predict  the  models  pattern.  However,  this  pattern  may 
differ  significantly  from  that  of  the  actual  corner  reflector. 

This  report  is  concerned  only  with  the  accuracies  of  moment  methods 
given  a thin-wire  model.  In  particular,  three  currently  available  and  often 
used  computer  programs  are  compared  by  their  application  to  simple  and 
complex  antenna  models.  Of  interest  are  not  highly  accurate  predictions  of 
current.  Instead,  "reasonable"  values  (within  'v  3 dB  of  converged  values) 
of  driving  point  impedance  and  radiation  patterns  (far-field)  are  sought. 

The  principal  variable  is  expansion  function  number  N>  and  a program  which 
generally  gives  satisfactory  results  for  least  N is  considered  superior. 

A number  of  thin-wire  moment  method  computer  program  comparisons  have 
been  conducted  previously  [Miller  1974,  Butler  1975,  Baldwin  1974,  Logan 
1973,  Rogers  1974].  However,  these  comparisons  either  (a)  considered  only 
scattering,  (b)  sought  high  accuracy  in  current  which  requires  many  expansion 
functions  per  wavelength  X,  (c)  considered  relatively  simple  shapes  such  as 
the  straight  wire,  "L"  shaped  wire  (2-wire  junction),  or  crossed  wires  (.4- 
wire  junction),  or  (d)  a combination  of  the  above. 

In  many  applications,  such  as  shipboard  antenna  performance  prediction, 
the  complexity  of  the  problem  often  precludes  high  accuracy.  However,  less 
accurate  predictions  are  often  satisfactory.  Also,  the  principal  parameters 
are  usually  radiation  patterns  and  driving  point  impedance  — not  current. 
Therefore,  the  computer  program  comparisons  presented  here  should  be  of 
interest  to  many  designers  or  analyzers  of  antenna  environments.  With 


regard  to  shipboard  applications,  a general  assessment  of  moment  methods  Is 
available  [Buchanan  et  al..  1977], 

In  Sections  2 and  3 a review  of  thln-wlre  moment  method  modeling  Is 
presented.  In  Section  4 three  computer  programs,  differing  primarily  in 
choice  of  expansion  and  testing  functions,  are  discussed.  The  means  for 
comparing  these  programs  is  detailed  in  Section  5 7 Finally,  in  Section  6 
the  programs  are  compared  by  presenting  driving  point  impedance  and  radia- 
tion pattern  computations.  These  results  are  obtained  by  considering  four 
different  wire  models  — two  simple  and  two  complex  — each  with  increasing 
numbers  of  expansion  functions.  Computer  run  time  is  also  considered. 


2.  THIN  WIRE  ANTENNA  THEORY 


The  most  popular  formulations  for  solving  thin-wire  antenna  problems 
are  of  the  E-field  type.  The  wire  surface  current  J is  sought  such  that 
the  E-field  it  radiates  satisfies  the  condition 

f 0 on  perfect  conducting  surface  of  wire 


-tan 


(1) 


^ applied  field  at  location  of  excitation  (gap)  on  wire 


where  I signifies  "component  tangential  to  wire  surface."  Two  simplifi- 
tan 

cations  are  usually  invoked  prior  to  implementing  these  formulations.  First, 
the  thin-wire  approximation  is  made.  With  reference  to  Figure  1 this  ap- 
proximation assumes  (a)  there  is  no  circumferentially  directed  component  of 
J,  (b)  J is  replaced  by  a filiament  of  current  I located  along  the  wire 
axis  C' , (c)  equation  (1)  is  satisfied  along  a line  C parallel  to  the  wire 
axis  and  offset  a wire  radius  a from  the  axis,  and  (d)  equation  Cl)  is 
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applied  only  to  the  axial  component  of  E.  Points  along  C'  are  defined  by 
the  length  variable  V and  those  along  C by  £. 

The  second  simplification  often  used  is  that  the  ith  applied  field  is 

A 

specified  as  a finite  voltage  over  a small  gap  (delta  gap  excitation). 
Thus  (1)  becomes 

L(I,Jt)  - l V.6GI  - £.)  U) 

i 1 1 

where 

L(M)  - -E'i  C3) 

?< 

A 

the  are  the  locations  of  applied  voltage  excitations  referenced  as 

A 

shown  in  Figure  2 and  l denotes  a unit  vector  at  l oriented  parallel  to  C. 

Most  frequency  domain,  moment  method,  thin-wire  computer  programs  are 
based  on  one  of  two  forms  for  the  operator  L. 

A.  Potential  Integrodif ferential  Equation 

Here  L has  the  form  [Harrington  1968  (Chapter  4)) 


L(I,£)  - jwy 


*-**I(*')GU,£')dit'  - j 


diq*) 


dt. 


r*-G(i,£«)d*«  (4) 


where 


G(M')  - 


;-jkR 

4ttR 


(5) 


and  R is  the  distance  between  a "source"  point  at  V (&'  is  corresponding 
unit  vector)  on  the  wire  axis  C'  and  a "field"  point  at  H on  the  wire  surface 
path  C.  The  remaining  variables  have  their  usual  meanings. 

fl.  Pocklington1 g Equation 

Here  L has  the  form  [Miller  1974] 
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L(I,£) 


/V  A 


(6) 


— f : 

jWE  J 

C' 


3£' 
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A major  difference  between  these  two  formulations  is  that  (A)  involves 
derivatives  of  the  unknown  current  I and  (6)  involves  derivatives  of  the 
Green's  function  G.  For  a finite  wire  radius  (A)  and  (6)  can  be  derived 
from  one  another.  Thus  exact  solutions,  if  they  exist,  are  independent  of 
the  formulation  used  to  obtain  them.  However,  when  approximate  solutions 
are  sought,  as  via  the  moment  method,  the  particular  formulation  used  may 
have  a significant  effect  on  convergence,  computational  ease,  etc. 

A third  form  for  L results  from  an  extension  of  Hallen's  integral  equa- 
tion [Butler  1972].  Hallen's  equation  is  devoid  of  derivatives,  thus  elimi- 
nating associated  computational  difficulties.  This  form,  however,  is  compli- 
cated by  the  introduction  of  unknown  constants  that  need  solving  along  with 
the  current.  Also  Butler's  form,  although  permitting  wire  juntions,  is 
restricted  to  straight  wire  configurations.  Other  modifications  of  Hallen's 
form  applicable  to  curved  wires  [Mei  1965]  and  curved  wires  with  junctions 
[Taylor  1969]  suffer  from  additional  operator  complexities.  For  perhaps 
these  reasons  t lie  authors  are  not  aware  of  a user-oriented,  arbitrary-wire- 
geometry  computer  program  based  on  a Hallen's  type  integral  equation.  This 
is  not  a severe  omission,  however,  for  the  following  reason  . Equation  C2) 
must  be  approximated  by  a linear  system  of  equations  before  it  can  be  solved 
for  arbitrary  geometries.  This  approximation,  the  crux  of  the  moment  method 
which  will  be  discussed  shortly,  involves  integrating  both  sides  of  C2)  after 
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multiplication  by  a testing  (or  weighting)  function.  There  is  a wide 
class  of  testing  functions  to  choose  from,  and  it  is  shown  by  Wilton  and 
Butler  [Wilton  1976]  that  a particular  weighting  of  Poeklington's  Equation 
is  equivalent  to  another  weighting  of  Hallen's  Equation. 


3.  METHOD  OF  MOMENTS 


The  method  of  moments  (or  moment  method)  is  a unifying  concept  for  re- 
ducing a linear,  inhomogeneous  integrodif f erential  equation  such  as  (2)  to 
a set  of  simultaneous,  linear,  algebraic  equations.  In  matrix  form  the 


latter  becomes 


[ZJI  = V 


where  [Z]  is  a known  square  matrix  and  V and  I are  known  and  unknown  column 
vectors  respectively.  An  in-depth  treatment  of  the  moment  method  is  given 
by  Harrington  [Harrington,  1968],  Once  suitable  [Z]  and  V are  known  (7)  can 
be  solved  for  I.  The  elements  of  I are  the  coefficients  in  an  expansion  of 
the  solution  to  the  original  equation.  Only  V contains  information  on  the 
sources  of  the  original  equation.  Thus  [Z]  need  be  computed  only  once  to 
solve  a time-harmonic  proolem  of  a particular  geometry  and  frequency  but 
where  the  source  distribution  varies. 

Equation  (7)  is  of  a form  often  appearing  in  network  theory  as  relating 
voltages  and  currents  of  an  N-port  network.  Thus  [Z]  is  often  referred  to 
as  the  "generalized  impedance"  matrix  and  I and  V the  "generalized  current" 
and  "generalized  voltage"  vectors  respectively. 

The  process  of  reducing  (2)  to  (7)  begins  with  the  selection  of  a set 
of  N Independent  functions  {f^}  defined  on  C*  with  which  to  expand  I.  Thus 
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(8) 


l. 


i 


N 

i(fc')  -Iif  on 

n-1  n n 

where  the  coefficients  1^  are  unknown.  A set  of  N independent  testing 

functions  {g  } defined  on  C are  then  chosen.  Finally  a suitable  inner 
m 

product  [Harrington  1968]  usually  of  the  form 


<<J>1 » 4>2> 


J ♦1CA)4>2C*)dJl 
c 


(9) 


is  chosen  and  the  inner  product  of  each  side  of  (2)  is  taken  after  multi- 
plication of  both  sides  with  each  g . The  resulting  N equations  in  N un- 

m 

knowns  (the  I ) are  expressed  in  matrix  form  by  (7)  where  the  nth  element 
n 

of  I is  I , the  rath  element  of  V is 
n 

Vm  [ V.«l>  «°> 


and  the  m,nth  element  of  [Z]  is 

Z 

mn 


j gma)L(fn,*)dJi 

c 


(11) 


3. 1 Basis  Functions 

The  Expansion"  or  "basis"  functions  fR>  with  which  I is  expanded,  must 

be  independent.  Also  the  f^  should  be  chosen  to  well  approximate  a basis 

for  the  domain  of  L.  One  expects  that  f individually  exhibiting 

n 

characteristics  of  1 will  result  in  fewer  of  them  and  consequently 
less  effort  to  solve  (7).  Thus  the  f^  are  usually  chosen  to  force  I to 
satisfy  current  continuity  at  each  point  along  C*.  This  implies  that  the  f 

n 

should  be  continuous  along  a wire  excluding  a junction  point  formed  by 
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- 


0 at  a wire  end.  Also  at  a multi-wire 


three  or  more  wires  and  that  f 

n 

junction  the  from  each  of  the  wires  are  often  distributed  such  that  the 
total  junction  current  satisfies  continuity. 

Occasionally  the  f are  not  within  the  domain  of  L.  For  example  pulse 
functions  (constant  over  a small  region  of  C’  and  zero  elsewhere)  are  not 
within  the  domain  of  L as  given  by  (4).  However,  such  simple  basis  functions 
can  still  be  used  if  L is  "approximated",  for  example,  with  a finite- 
difference  evaluation  of  the  derivatives.  The  computer  program  WRSMOM,  which 
is  discussed  later,  is  based  on  pulse  expansion  functions  with  a finite- 
difference  approximation  to  L.  Further  discussions  concerning  "approximate 
operators"  and  other  techniques  such  as  "extended  domains",  etc.  are  given 
in  [Harrington  1968]. 

The  f^  are  generally  composed  of  one  of  two  function  types  — entire 
functions  or  subsectional  functions.  Subsectional  functions  are  each  de- 
fined other  than  zero  only  over  a small  region  of  C'.  An  example  is  the 
pulse  functions  discussed  above.  Entire  functions  are  not  so  constrained. 
Examples  of  entire  function  bases  are  Maclaurln  Series  polynomials,  Fourier 
trigonometric  functions  and  Legendre  polynomials.  These  three  bases  are 

comprised  of  independent  functions.  In  addition,  the  latter  two  are  com- 
prised of  functions  orthogonal  with  respect  to  (9).  However,  neither  basis 

necessarily  orghogonalizes  L (<4>^,L<^2>  = 0 if  T*  which  would  result  in  di- 
agonal [Z],  A set  of  functions  which  does  diagonalize  [Z]  would,  of  course, 
form  an  excellent  basis  but,  for  arbitrary  geometries  such  a set  of  functions 
is  costly  to  compute.  (The  class  of  rotationally  symmetric  bodies  is  an  excep- 
tion in  that  basis  functions  with  sinusoidal  circumferential  variation  are 


known  to  "block"  orthogonalize  L [Harrington  1969,  Mautz  1977].) 

Experience  indicates  that  entire  function  expansions  tend  to  converge 
faster  than  subsectional  function  expansions  [Thiele  1973].  However,  the 
latter  tends  to  result  in  "better  conditioned"  [Z]  matrices  in  that  the  ratio 
of  diagonal  to  off-diagonal  [Z]  elements  has  generally  greater  magnitude.  A 
procedure  for  solving  (7)  when  [Z]  is  better  conditioned  is  generally  less 
sensitive  to  round-off  errors. 

The  computation  of  [Z]  for  an  entire  function  expansion  is  usually 

considerably  more  involved  than  for  a subsectional  function  expansion.  Thus 

for  complex  wire  geometry  problems  a subsectional  function  expansion  is  more 

inviting  than  an  entire  function  expansion.  Hence,  the  vast  majority  of 

available  moment  method  computer  programs  for  modeling  arbitrary  collections 

of  bent  wires  are  based  on  subsectional  function  expansions.  Three  examples 

of  subsectional  functions  are  shown  in  Figure  3.  The  f in  3(a)  are  pulses 

n 


defined  by 


f (£')  = 
n 


1 £ -A  /2  < V < £ +A  /2 

n n n n 


0 Otherwise 


where  A is  the  pulse  width  defined  in  the  figure.  The  f in  3(b)  are  over- 
n n 


lapping  triangles  given  by 


i (£'-£  , )/A  £<£'<£ 

) n-1  n n-1  — — n 

. <Wl'K 


Otherwise 


where  An  and  are  defined  in  the  figure.  The  fR  in  3(c)  are  overlapping 


sinusoids  given  by 
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fn-l 


fn 


fn  + l 


(b) 


Figure  3.  Subsectional  functions:  (a)  pulse 

(piece-wise  constant),  (b)  triangular 
(piece-wise  linear),  and  (c)  piece-wise 
sinusoidal. 
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Geometry  description  and  [Z]  computation  in  a computer  program  appli- 
cable to  arbitrary  configurations  of  bent  wires  is  facilitated  usually  by 
approximating  the  paths  C'  (and  C)  with  straight  line  segments.  A typical 
segmentation  is  shown  in  Figure  4.  For  pulse  expansion  functions  the  seg- 
ments usually  coincide  with  the  Afl  of  Figure  3.  For  triangle  and  piece-wise 
sinusoid  expansion  functions  the  segments  generally  coincide  with  the  A^ 

(and  ^n). 

As  previously  mentioned,  it  is  desirable  to  incorporate  within  the  f 
known  constraints  on  I.  Continuity  of  1 is  one  such  constraint.  The  pulse 
functions  of  Figure  3(a)  do  not  enforce  continuity  of  I.  However,  they  are 
generally  used  in  conjunction  with  approximate  operators  that  afford  addi- 
tional smoothing.  The  triangle  and  sinusoid  (piece-wise)  functions  of 
Figure  l(b,c),  on  the  other  hand,  do  enforce  continuity  of  current  at  each 
point  along  (.'.  These  two  bases  also  maintain  continuity  of  current  at 
mut 1 ip le-wi re  junctions  if  a simple  rule  is  followed  in  modeling  the  wires. 
To  illustrate,  consider  the  j-wire  junction  shown  in  Figure  5(a).  The 
modeling  begins  by  choosing  any  one  of  the  wires  as  terminated  at  the  junc- 
tion (Figure  3(b)).  Then  the  remaining  wires  are  added  in  succession  such 
that  each  overlaps  any  previously  placed  wire  (Figure  5(c,d)).  This  over- 
lapping is  accomplished  by  aligning  the  junction  point  with  the  peak  of  the 
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end  triangle  (or  sinusoid)  function  on  the  wire  being  added.  The  justifica- 
tion for  this  method  is  given  in  [Chao  1970].  It  is  also  presented  in  terms 
of  "independent  loop  currents"  in  [Richmond  (no  date)]. 

3.2  Testing  Functions 

The  set  of  testing  functions  [g  } must  also  be  independent.  Other  re- 

m 

strictions  on  the  g are  not  obvious  from  their  explicit  use  as  indicated  by 
m 

(10)  and  (11).  However,  one  notes  that  the  moment  method  is  equivalent  to 

g 

the  Rayleigh-Ritz  procedure  if  the  gm  ire  in  the  domain  of  L the  adjoint 
operator  of  L [Harrington  1968]  defined  by 

<4>1,L4>2>  - <La4>2,4>1> 

This  equivalence  is  noteworthy  since  the  Rayleigh-Ritz  procedure  achieves  an 

approximate  solution  to  a variational  expression  for  a scalar  quantity  p of 

interest  (impedance,  far-field,  etc.),  i.e.  small  errors  in  I when  used  to 

compute  a variational  expression  for  p result  in  proportionately  less  error 

in  p.  For  the  "symmetric"  inner  product  defined  by  (9)  the  operators  of 

(4)  and  (6)  are  self-adjoint,  i.e.,  L = L.  Thus  it  is  desirable  to  choose 

the  g from  the  domain  of  L. 
m 

Two  types  of  {g^}  are  often  found  in  existing  moment  method  computer 

programs.  One  type  sets  {gm)  = {f^}.  This  is  often  called  Galerkin's 

method.  One  advantage  of  Galerkin's  method  is  that  [Z]  is  then  symmetric. 

This  decreases  the  effort  required  for  computing  [Z]  and  for  solving  (7). 

It  also  saves  computer  storage  since  only  the  upper  or  lower  triangle  of  [Z] 

is  needed.  A disadvantage  is  that  the  g may  complicate  the  evaluation  of 

m 

(11) .  Examples  of  programs  based  on  Galerkin's  method  are  the  Syracuse  and 


15 


Ohio  State  programs  (subsectional  triangles  and  subsectional  sinusoids  re- 
spectively) discussed  later. 


The  other  frequently  chosen  {g  } is  the  set  of  impulses  positioned  at 

m 

the  centers  of  the  f . This  choice  of  g reduces  the  integral  computation 

n m 

in  (11)  to  a trivial  operation.  The  choice  of  impulse  functions  for  testing 

is  also  appropriately  called  "co-location"  or  "point-matching".  Of  course, 

these  g are  not  within  the  domain  of  L.  However,  they  are  often  used  in 
m 

conjunction  with  a finite-difference  approximation  to  L which,  as  discussed 

earlier,  provides  a compensating  smoothing  effect.  In  fact,  it  has  been 

shown  that  a finite-difference  approximation  to  Pocklington's  equation  is 

closely  akin  to  testing  with  triangles  or  sinusoids  [Wilton  1976],  The 

computer  program  WRSMOM  discussed  later  uses  impulse  testing. 

Equation  (10)  takes  a different  form  when  g is  an  impulse  function. 

m 

In  this  case  the  ith  antenna  excitation  is  generally  treated  as  a specified 
voltage  VA  over  a small  but  finite  gap.  The  gap  size  is  equal  to  the  seg- 
ment width  A of  the  n th  pulse  expansion  (basis)  function  occurring  at 
ni 

the  ith  excited  port.  Then  (2)  becomes 


L(I-*'>  -I’iF 


The  testing  functions  are  weighted  impulses  of  the  form 


g - 6(*  - i)Am 
m mm 


(15) 


(16) 


Thus  under  impulse  testing  and  pulse  expansion  (10)  and  (11)  are  replaced  by 


if  the  ith  excitation  occurs 
at  mth  expansion  function  (m*n 

Otherwise 


i> 

(17) 


- L(f  ,Z  )A 
mn  ' n’  m m 


(18) 
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* 9m 


a)  FOUR-PULSE  APPROXIMATION  TO  fn(-E') 
AND  FOUR-IMPULSE  APPROXIMATION 
TO  gm(  -E  ) 


b) 


FOUR-PULSE  APPROXIMATION  TO  df"(f  } 

dz 

AND  FOUR-IMPULSE  APPROXIMATION 


TO 


dgmd) 
d l 


Figure  6.  Syracuse  program  (SYR)  treatment  of 
(a)  expansion  and  testing  functions, 
and  (b)  their  derivatives. 
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oriented  programs  of  general  applicability  based  on  the  pulse-expansion  im- 
pulse-testing method  originally  proposed  by  Harrington  [Harrington  1967]. 

The  expansion  functions  are  defined  by  (12)  and  the  testing  functions  by  (16). 
The  operator  L is  given  by  (4).  However,  as  discussed  in  the  previous  sec- 
tion, a finite-difference  approximation  to  L is  employed. 

The  expansion  and  testing  functions  in  WRSMOM  are  perhaps  the  simplest 
feasible.  Thus  [Z]  is  relatively  simple  to  compute.  However,  [Z]  is 
generally  not  symmetric. 

The  principal  drawback  to  WRSMOM  is  that  the  pulse  expansion  functions 
do  not  constrain  current  continuity  at  not  only  multiwire  junctions  but  also 
at  segment-to-segment  interfaces  along  a single  wire.  However,  as  previously 
mentioned,  use  of  the  approximate  form  of  L tends  to  "smooth"  these  discon- 
tinuities. This  has  been  demonstrated  [Mittra  1975]  if  the  segment  sizes  are 

uniform.  But  uniform  A is  a severe  restriction  if  large,  complex  structures 

n 

are  involved.  Thus  A variations  in  excess  of  2:1  have  been  permitted  in 

n 

some  of  the  examples  discussed  in  Section  6. 

A WRSMOM  program  description  is  available  [Warren  1974], 

5.  PROGRAM  COMPARISON  CONSIDERATIONS 
Current  itself  is  usually  not  of  ultimate  interest  in  the  prediction  or 
analysis  of  antenna  performance.  Radiation  pattern  and  driving  point  imped- 
ance are  usually  the  desired  parameters.  Once  the  current  has  been  found, 

*Impulse  expansion  functions  (Hertzian  dipoles),  of  course,  are  even  simpler 
but  they  have  shown  promise  only  in  "hybrid"  programs.  These  programs  use 
impulse  expansion  functions  only  for  computing  [Z]  elements  corresponding  to 
distant  current  elements  [Stewart  1974  (Section  3C) ] . 
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however,  e.g.,  by  the  moment  method,  patterns  and  impedances  are  easily  ob- 
tained. The  three  programs  described  in  the  last  section  each  compute 
patterns  and  impedances.  Therefore,  a meaningful  comparison  between  these 
programs  is  best  made  by  assessing  their  relative  ease  in  predicting  patterns 
and  impedances  to  within  a desired  accuracy  for  a variety  of  structures. 

The  accuracy  of  a method  for  predicting  the  EM  characteristics  of  a 
complicated  structure  is  difficult  to  assess.  This  is  because  of  the 
scarcity  of  either  exact  solutions  or  reliable  experimental  results  for 
complicated  structures.  Therefore,  "correct"  solutions  are  defined  here  as 
those  values  that  all  program  results  converge  to  that  are  reasonable.  Pro- 
gram results  converge  if  the  addition  of  expansion  functions  does  not  ap- 
preciably change  results. 

Once  converged  parameter  values  (patterns  and  impedances)  are  obtained, 
a program  is  evaluated  by  its  application  with  fewer  and  fewer  expansion 
functions.  A minimum  number  of  expansion  functions  is  determined  when  signi- 
ficant 3uB)  deviations  from  the  converged  values  result.  The  program 
exhibiting  smallest  minimum  is  considered  superior  for  that  structure  type. 

Of  course,  consideration  must  also  be  given  to  other  variables  contributing 
to  computer  run  time  and  storage  requirements.  For  example  if  OS  requires 
more  expansion  functions  than  WRSMOM,  it  does  not  necessarily  mean  that  the 
run  time  for  computing  [Z]  ("matrix  fill"  time)  is  less  for  WRSMOM  than  for 
OS.  Nor  does  it  mean  that  the  run  time  for  solving  (7)  or  even  that  the  [Z] 
storage  requirements  are  less  for  WRSMOM  than  for  OS.  This  is  because  the 
symmetric  matrix  characteristic  of  OS  gives  it  a significant  advantage  at 
the  outset.  Thus,  although  expansion  function  number  is  the  principal 
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variable  for  program  comparisons  of  the  following  section,  the  overall  evalu- 
ation considers  actual  timing  and  storage. 

Driving  point  impedance  computations  are  considered  within  desired  ac- 
curacy if  both  real  and  imaginary  parts  differ  from  converged  values  by  at 
most  a factor  of  ^ 2 (or  ^1/2). 

Radiation  gain  patterns  are  considered  acceptable  if  they  differ  from 
converged  patterns  by  at  most  ^ 3dB  at  all  pattern  points. 

In  all  cases  the  driving  point  impedance  Z is  computed  as  the  ratio 

A 

of  the  applied  voltage  V,  to  the  resulting  coefficient  I of  current  basis 

1 nl 

function  at  the  location  of  the  source.  The  radiation  gain  patterns  are  ob- 
tained by  dividing  4tt  times  the  far-field  intensity  (for  a desired  polariza- 
tion, 0 or  0,  of  E-field)  by  the  average  input  power  P.  Here  P is  determined 
from 

P - IVJ2  Real{l/Zin) 

A 

Note  that  is  an  RMS  value. 

6.  PROGRAM  COMPARISONS 

Results  of  applying  the  three  computer  programs  OS,  SYR,  and  WRSMOM  to 
four  thin-wire  structures  are  discussed  here.  The  first  two  structures  con- 
sidered are  relatively  simple:  straight,  resonant  wire  and  three-wire  junc- 
tion. The  second  two  are  complex:  whip-excited  twin  fan  and  self-excited 
twin  fan.  These  structures  are  defined  below. 

As  mentioned  in  Section  1 program  comparisons  with  application  to  simple 
structures  such  as  a straight  wire  and  single  three-  or  four-wire  junctions 
have  previously  been  reported  on.  The  emphasis,  however,  was  usually  on 
obtaining  a high  degree  of  accuracy  in  current.  Here,  on  the  other  hand. 
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the  goal  Is  to  compare  the  computational  efforts  needed  to  obtain  reasonable 
far-field  and  Impedance  values. 


The  coordinate  system  referred  to  throughout  this  section  is  shown  in 
Figure  7.  Also  the  variable  for  "number  of  expansion  functions"  is  denoted 
N. 

Except  where  otherwise  mentioned,  OS  was  run  with  the  4-interval 
Simpson's  rule  integration  discussed  in  Section  4.  The  exception  pertains 
to  the  self-excited  twin  fan  example. 

All  examples  have  uniform  wire  radii  since  OS  is  not  otherwise  appli- 
cable. Of  course,  this  could  be  a serious  program  constraint.  However,  as 
previously  mentioned  in  Section  4,  the  many  simplifying  assumptions  usually 
required  in  reducing  a complex  real-world  structure  to  a tractable  wire 
model  are  likely  to  dwarf  that  of  uniform  wire  radius. 

The  distribution  of  expansion  functions  throughout  a wire  model  was 
maintained  approximately  uniform  whenever  practical.  A noted  exception  is 
the  self-excited  twin  fan  example  where  the  excited  wires  are  many  times 
smaller  than  any  of  the  other  wires.  This  extreme  wire-length  variation 
provided  a severe  test  for  the  three  programs  and  thus  served  a purpose  of 
this  report.  However,  such  difficulties  can  often  be  avoided  by  proper 
modeling.  An  example  of  a better  model  for  the  self-excited  twin  fan  is 
discussed  later.  The  self-excited  twin  fan  example  also  violates  the  usual 
thin-wire  criterion  that  the  length  L of  each  wire  be  many  times  the  wire 
radius  a.  The  excited  wires  have  L/2a  - 6. 

The  whip-excited  twin  fan  model  has  junctions  with  acute  angles  between 
wires  as  small  as  ^lS®.  This  also  provides  a severe  test  for  the  thin-wire 
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programs.  Acute  angles  less  than  30°  are  purported  to  be  avoided  [Richmond 
(no  date) ] . 

Many  of  the  results  were  obtained  with  a minimum  of  expansion  functions. 
Thus  unusually  large  segment  lengths  A occurred.  This  necessitated  bypassing 
a check  incorporated  in  OS  for  "too-long"  A's. 

The  A/a  ratios  were  maintained  significantly  greater  than  unity.  This 
is  important  since  moment  method  solutions,  with  sinusoidal  expansion  func- 
tions, of  problems  modeled  with  and  without  the  thin-wire  approximation  have 
been  demonstrated  to  deviate  significantly  for  A/a  = 1 [Imbriale  1973],  In 
fact  a diverging  nature,  not  attributable  to  round-off  error,  of  the  solu- 
tion with  decreasing  A/a  was  observed.  In  view  of  this  and  if  it  can  be 
assumed  that  a subsectional  basis  is  complete,  then  it  follows  that  no  solu- 
tion exists  for  models  incorporating  the  thin-wire  approximation.  Of  course, 
approximate  solutions  such  as  those  obtained  with  A/a  > 1 apparently  do  re- 
semble the  solution  to  actual  wire  problems.  Thus  solutions  for  thin-wire 
approximated  models  appear  to  asymptotically  approach  solutions  for  the 
actual  wire  models  in  that  A/a  must  be  neither  too  large  nor  too  small, 

6.1  Simple  Structures 
A.  Resonant  straight  wire 

This  center-excited  structure  is  shown  in  Figure  8.  The  frequency  is 
300  MHz.  Table  I gives  the  driving  point  impedances  computed  by  the  three 
programs  for  different  numbers  of  expansion  functions  (N’s).  Equal  length 
segments  were  maintained.  For  convenience,  only  odd  N's  were  chosen  since 
the  wire  is  centrally  excited.  Reasonable  agreement  between  the  three 
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WIRE  RADIUS  = I0"4  m 
EXCITATION  = 1.0  V 
FREQUENCY  * 300  MHz 


Figure  8.  Resonant  straight  wire. 


Dashed  line  underlines  impedances  which  first  fall  within 
convergence  criteria  discussed  in  Section  5. 
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programs  Is  apparent  with  N ■ 21  and  the  corresponding  values  were  assumed 
correct.  The  dashed  line  underlines  the  Impedances  which  first  fall  within 
convergence  criteria  discussed  In  Section  5.  Thus  OS  demonstrates  better 
convergence  here.  The  OS  program  does  not  permit  an  N - 1 run.  However, 
it  Is  expected  that  the  single  expansion  function  case  for  OS  would  also 
give  acceptable  Impedance  values  since  it  is  simply  the  usual  sinusoidal- 
current  assumption  [Jordan  1950  (pg  543)].  On  the  other  hand,  it  is  possible 
that  SYR  and  WRSMOM  runs  with  N = 5 would  also  be  acceptable.  These  were 
not  obtained  at  the  time  of  data  taking. 

In  Figures  9-11,  X-Y  plane  radiation  gain  patterns  are  given  for  the 
component  of  E-field.  The  three  programs  display  excellent  results  for  the 
minimum  number  of  expansion  functions  that  can  be  applied. 

B.  Three-wire  1 unction 

This  antenna  is  shown  in  Figure  12.  The  total  wire  length  is  'U).7A. 
Approximately  equal  segment  lengths  were  maintained  throughout.  In  Table  II 
driving  point  Impedances  computed  by  the  three  programs  again  exhibit 
fastest  convergence  for  OS.  The  minimum  number  of  expansion  functions  N = 3 
appears  adequate  for  OS.  This  number  is  arrived  at  by  centering  one  expan- 
sion function  at  the  excitation  and  providing  two  more  to  satisfy  current 
continuity  at  the  junction  as  discussed  in  Sections  3 and  4. 

In  Figures  13-15  the  X-Y  plane  radiation  gain  patterns  for  the  0 compo- 
nent of  E-field  are  given.  All  three  programs  provide  acceptable  (within 
3 dB  of  converged  pattern,  as  discussed  in  the  previous  section)  results 
with  the  minimum  number  of  expansion  functions  N « 3.  However,  the 
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Radiation  gain  pattern  (if -component  of  E-field)  convergence  for  OS  applied  to  resonant  wire 


3.0  45.0  90.0  135.0  180.0  225.0  270.0  315.0  360.0 

PHI  (DEGRrES) 

Figure  11.  Radiation  gain  pattern  (1)1 -component  of  E-field)  convergence  for  WRSMOM  applied  to 


Figure  12.  Three-wire  Junction. 
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TABLE  II 


Driving  Point  Impedance  for  Three  Wire  Junction 


No.  of  Expansion 
Functions  (N) 


OS 


15* 


82. 1+j  7. 2 


3 

6 

** 

9 

II*** 


75 . 0+j  3. 7 
78. 7+j4. 8 
80. 9+j  5. 8 
81 . 2+ j6. 7 


Program 

SYR 

WRSMOM 

79. 0+j  3 . 0 

75. 0+j 13.0 

46. 3- j51 . 3 

69.4+j47.4 

66. 6-j 11 . 2 

87 , 7+j24. 3 

75 . 5- j 2 . 0 

71.4+jl8.4 

76. 6+ j 1 . 3 

89.3+j 12.6 

Dashed  line  underlines  impedances  which  first  fall  within  convergence 
criteria  discussed  in  Section  5. 


* 13  for  WRSMOM 

**  8 for  WRSMOM 

***  10  for  WRSMOM 
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convergence  fox  OS  (Figure  13)  is  exceptionally  rapid. 

In  Figures  16-18  the  4>  component  of  E-field  radiation  gain  patterns 
are  plotted.  Again  the  convergence  for  OS  is  exceptional  (Figure  16),  and 
the  minimum  N for  WRSMOM  is  adequate  (Figure  18).  However,  SYR  appears  to 
require  additional  expansion  functions  to  achieve  satisfactory  results 
(Figure  17). 

It  is  interesting  that  WRSMOM  demonstrates  at  least  as  good  convergence 
as  SYR  for  this  problem.  Logan  [Loganl973]  also  applied  WRSMOM  and  SYR  to 
a three-wire  junction  problem  (albeit  a configuration  different  from  that 
of  Figure  12).  He  found  WRSMOM  to  be  significantly  less  accurate  than  SYR 
which  he  attributed  to  the  presence  of  the  junction.  However,  a discussion 
with  WRSMOM's  author  [Warren  1975]  suggests  that  the  presence  of  a large 
adjacent  subsectional  variation  (>  2:1)  is  the  real  culprit  in  Logan's 
comparison.  In  Section  4 it  was  pointed  out  that  the  pulse  expansion  func- 
tions used  in  WRSMOM  do  not  constrain  current  continuity  even  on  single 
wires.  Thus  the  junction  presence  alone  should  not  cause  large  errors. 

6.2  Complex  St»  ctures 
A.  Whip-excited  twin  fan 

This  structure  is  shown  in  Figure  19.  The  whip  antenna  is  5m  from  the 
simplified  twin-fan  model.  A ground  plane  is  assumed  and  thus  the  presence 
of  the  Images.  However,  for  generality  the  computer  programs  do  not  in- 
corporate the  simplifications  made  possible  by  this  symmetry.  The  excita- 
tion port  (at  the  6-wire  junction)  of  the  twin-fan  is  assumed  shorted. 

A total  of  18  wires  comprises  this  model  which  includes  two  6-wire 
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Figure  16.  Radiation  gain  pattern  -component  of  E-field)  convergence  for  OS  applied  to 
three-wire  junction  antenna. 
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three-wire  junction  antenna. 


junctions.  Its  complexity  necessitated  convenient  geometry  data  Input 
checks.  Thus  each  of  the  three  programs  was  modified  to  generate  four  plots 
of  the  geometry  being  modeled  after  Input  data  Is  read.  The  four  plots  are 
the  views  shown  in  Figure  19.  This  provided  an  alert  to  data  input  error. 

In  Table  III  the  driving  point  impedances  are  compared.  Again  OS  ex- 
hibits superior  convergence.  However,  SYR  and  WRSMOM  converge  to  reasonable 
values  almost  as  rapidly.  The  minimum  number  of  expansion  functions  (N « 18) 
for  OS  and  SYR  were  arrived  at  by  requiring  one  for  the  excitation  and  ad- 
ditional ones  to  meet  the  junction  current  continuity  constraint  (at  the 
2-  as  well  as  6- wire  junctions)  as  discussed  in  Sections  3 and  4.  The 
N - 18  case  for  WRSMOM  simply  assigned  an  expansion  function  to  each  wire. 
This  resulted  in  severe  adjacent  segment  length  variations. 

The  manner  in  which  expansion  functions  for  the  N » 80  case  for  all 
three  programs  was  distributed  assumed  - '''•10  expansion  functions  per  wave- 
length throughout  the  model.  However,  the  adjacent  segment  size  variation 
was  still  in  excess  of  2:1.  As  mentioned  in  Section  6.1  this  variation 


could  lead  to  errors  in  WRSMOM  generated  results.  However,  it  is  Important 
to  permit  such  variations  if  large,  complex  structures  are  to  be  modeled. 

Radiation  gain  patterns  for  the  9 component  of  E-field  in  the  X-Y  plane 
are  given  in  Figures  20-22  for  "tight"  convergence  (large  N).  Results  from 
all  three  programs  converge  to  approximately  the  same  patterns.  Note  Y-axis 
symmetry  (about  <p  “ 90°  and  about  <p  = 270°)  as  expected. 

A "loose"  convergence  is  shown  in  Figures  23-25.  OS  and  WRSMOM  dis- 
play satisfactory  convergence  with  N ■ 21  but  SYR  requires  N ■ 28.  It  is 
Interesting  that  WRSMOM  shows  a little  better  convergence  than  even  OS  here. 
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TABLE  III 

Driving  Point  Impedance  for  Whip  Excited  Fan 


No.  of  Expansion 
Functions  (N) 

OS 

Program 

SYR 

WRSMOM 

80 

48.2-jl22.0 

46 . 7- j 127 . 6 

18 

32. 1-j 125. 0 

24.8-j219. 0 

23.4-j42.1 

21 

45.9- j 122. 0 

33.5-2159.0 

42 . 4- j 1 14 . 8 

22 

48.2-j 122.0 

3i.9-jl57.0 

39.5-jll3.3 

24 

48.2-jl22.0 

30.4-jl55.0 

39. 6-j 113.3 

28 

48.2-jl22.0 

41 . 8-j 150. 0 

42. 1-j 113. 2 

52 

47. 7-j 123. 0 

45. 6-j 138. 0 

45 . 3- j 124 . 8 

68 

48. 1-j 123. 0 

46.7-jl32.0 

45. 7-j 127.4 

Dashed  line  underlines  impedances  which  first  fall  within  convergence 
criteria  discussed  in  Section  5. 


G.O  45. G 90. G 135. G 18C.G  225. G 270. G 315. G 36G.G 

Prir  (DEGREES) 

Figure  21.  Radiation  gain  pattern  (0-component  of  E-field)  convergence  (tight)  for  SYR  applied  to 
whip-excited  twin  fan. 
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Figure  25.  Radiation  gain  pattern  (9-component  of  E-field)  convergence  (loose)  for  WRSMOM  applied 
to  whip-excited  twin  fan. 


B.  Self-excited  twin  fan 


This  example  provides  the  severest  test  for  the  three  programs.  This 
is  readily  apparent  in  Figure  26  which  depicts  blow-ups  of  the  compli- 
cated feed  regions  (excitation  point  and  its  image).  The  feed  wires,  or 
connecting  arms,  separate  two  multiwire  junctions  (five-wire  and  three-wire) 
by  only  'V2.5xlO~^A.  The  model  is  otherwise  the  same  as  the  twin-fan  of 
Figure  19.  The  whip  is  omitted.  Thus  there  are  a total  of  19  wires  on 
this  model. 

"Tight"  convergences  of  radiation  gain  patterns  are  shown  in  Figures 
27-30.  The  OS  patterns  of  Figure  27  were  obtained  with  the  4-interyal 
Simpson's  Rule  integration  option.  These  patterns  were  repeated  with  the 
exponential  integral  option  with  results  given  in  Figure  30.  A comparison 
of  Figures  28  and  30  shows  that  patterns  generated  by  OS  with  the  exponential 
integral  option  and  SYR  converge  to  the  same  values.  Thus,  for  this  example, 
the  remainder  of  the  patterns  and  all  the  impedances  generated  by  OS  em- 
ployed the  exponential  integral  option.  The  patterns  using  WRSMOM  (Figure 

29)  are  unexpectedly  low  in  average  levels.  Thus  severe  errors  are  be- 
lieved to  occur  with  WRSMOM  applied  to  this  model.  Figures  31-33  compare 
OS,  SYR  and  WRSMOM  generated  patterns  when  the  feed  wires  (connecting  arms) 
are  increased  in  length.  The  three  programs  now  compare  reasonably  well  as 
expected  and  the  resulting  patterns  approximate  the  converged  patterns 
computed  by  OS  and  SYR  with  unlengthened  connecting  arms  (Figures  28  and 

30) .  Thus  the  latter  patterns  are  taken  to  be  "correct".  Driving  point 
impedances  were  also  computed  for  these  different  feed  wire  lengths.  These 
Impedances,  shown  in  Table  IV,  bear  out  the  above  conclusions. 
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Figure  27.  Radiation  gain  pattern 
self-excited  twin  fan. 
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TABLE  IV 

Fan  Impedance  for  Different  Excitation  Arm  Lengths 


Excitation 

Arm  Length  (x) 

No.  of 

Exp.  Fns.* 

OS 

Program 

SYR 

WRSMOM 

2.44xl0-3 

75 

49.4-j206.3 

52 - 2-j  211.7 

303.1-j2,915.0 

0.03 

75 

61.1-j207.4 

67. 5-j  215. 5 

44.5-j306.1 

0.06 

79 

69. 3-j 182 . 2 

77.0-jl95.6 

o7.2-j384.3 

0.12 

79 

1 

80. 4-j  93.2 

88.7-jl09.5 

50.7-jl05.7 

* Fot  WRSMOM  subtract  1 from  these  numbers. 


"Loosely"  converging  patterns  are  shown  in  Figures  34-36.  The  conver- 
gence for  OS  Is  a bit  superior  to  that  for  SYR  In  that  N = 27  provides  ac- 
ceptable results  for  OS  (Figure  34)  and  N = 29  for  SYR  (Figure  28)..  Although 
the  WRSMOM  convergence  Is  excellent  (Figure  36),  there  is,  as  discussed  above, 
considerable  doubt  as  to  its  accuracy. 

The  convergence  for  driving  point  impedance  is  demonstrated  in  Table  V. 
Both  OS  and  SYR  require  the  same  N to  achieve  acceptable  results. 

The  difficulty  with  WRSMOM  for  this  example  and  the  need  for  the  more 
time  consuming  exponential  integral  option  in  OS  can  probably  be  alleviated 
by  modeling  the  Junction  differently.  For  example  the  connecting  arms  be- 
tween junctions  can  be  removed  by  "passing  the  source  through  the  junctions" 
such  that  Kirchoff's  Voltage  Law  remains  satisfied  [Burke  1976].  Since 
two  junctions  are  involved,  the  equivalent  representation  is  not  unique. 

One  representation  is  shown  in  Figure  37.  Here  1/2  the  source  is  passed 
through  each  junction  and  the  short  connecting  arm  is  then  eliminated.  The 
driving  point  impedance  can  now  be  approximated  by  dividing  the  algebraic 
sum  of  the  currents  entering  the  junction  from  points  E and  F into  1.0V. 


6.3  Computer  Run  Time 

2 

The  matrix  fill  ( [ Z ] computation)  time  is  proportional  to  N and  the 

3 

simultaneous  equation  solution  time  is  proportional  to  N . Thus  the  compu- 


ter cpu  time  T in  seconds  for  solving  a moment  method  problem  is  given  by 

2 3 

T - AN  + BN 


for  reasonably  large  values  of  N,  e.g.  N > 40.  The  coefficients  A and  B 
for  the  three  programs  considered  are  given  below.  Values  for  the  Honeywell 
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Figure  34.  Radiation  gain  pattern  (■9-component  of  E-field)  convergence  (loose)  for  OS  applied  to 
self-excited  twin  fan.  (Exponential  integral  form). 
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0.0  45.0  90.0  135.0  180.0  225.0  270.0  315.0  360.0 
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Figure  35.  Radiation  gain  pattern  (8-component  of  E-field)  convergence  (loose)  for  SYR  applied  to 
self-excited  twin  fan. 


PHI  (DEGREES) 

Radiation  gain  pattern  (0 -component  of  E-field)  convergence  (loose)  for  WRSMOM  applied 
to  self-excited  twin  fan. 


TABLE  V 

Driving  Point  Impedance  for  Self  Excited  Fan 


No.  of  Expansion 

Functions  (N)*  OS 


75  49.4-j206.3 


20 

8.2-j297.1 

21 

6. 2-j  300.3 

23 

6.5-j300.3 

27 

48.6-j214.8 

28 

47.3-j216.4 

29 

48. 3-j217. 3 

51 

47.9-j209.4 

65 

49.4-j206. 5 

Program 

SYR 

52 . 2-j  211.7 

234. 1- j653.4 

84.4- j239.8 
79. 7-j 158. 5 

103.7:227^2 

114.2- j253.4 
108.8-j238. 0 

58.9-j222.5 

59.5- j219.2 


WRSMOM 


359. 8- j7, 591.0 

338. 1- j7, 682.0 
303. 3-j7, 743.0 
454. 6-j7, 664.0 

519^4:26,627^0 

403. 9- j4, 164.0 

243. 1- j3, 068.0 


Dashed  line  underlines  impedances  which  first  fall  within  convergence 
criteria  discussed  in  Section  5. 


*For  WRSMOM  subtract  1 from  these  numbers. 
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Figure  37.  Ce)  Excited  complex  junction. 

(b)  Equivalent  source  representation. 


(Multics)  6000  series  system,  the  computer  used  throughout  this  work,  and 
also  for  the  CDC  7600  [Miller  1974]  are  shown.  The  WRSMOM  values  for  the 
CDC  system  were  not  determined  from  WRSMOM  but  from  a Lawrence  Livermore 
Laboratory  program  written  to  accomplish  the  same  thing.  Also  the  OS  value 
of  A in  parentheses  corresponds  to  use  of  the  closed-form  exponential  inte- 
gral option  for  computing  [Z ] . The  other  values  assume  the  4-interval  Simp- 
son's Rule  option. 


Honeywell 

. (Multics)  6000 

CDC  7600 

OS 

A 

9.36 

X 

10~3  (4.54  x 10~2) 

3.5  x 10"4 

B 

1.66 

X 

m 

l 

o 

rH 

1.07  x 10" 

SYR 

A 

5.81 

X 

10-3 

-4 

2.6  x 10 

B 

9.18 

X 

10-5 

5.0  x 10‘6 

WRSMOM 

A 

9.09 

X 

io-3 

3.0  x 10~4 

B 

9.18 

X 

10-5 

5.3  x 10~6 

OS  employs  the  square  root  method  [Faddeey  1963  (pp. 144-147) ] for 
solving  (7).  This  is  considerably  faster  than  the  general  matrix  inversion 
method  Incorporated  in  both  SYR  and  WRSMOM.  However,  one  can  easily  fit  to 
either  SYR  or  WRSMOM  a Gaussian  elimination  scheme  which  simply  solves  (7). 
This  would  reduce  their  B values  by  about  1/3.  Also  SYR  is  approximately 
symmetric.  Thus  incorporation  of  the  square  root  method  in  SYR  might  be 
feasible.  This  would  further  reduce  B (and  alao  lower  A)  for  SYR. 

For  equal  N the  matrix  fill  time  for  SYR  is  smaller  than  for  OS.  This 
is  especially  true  if  the  latter  is  run  with  the  accurate  but  slow  exponen- 
tial integral  option.  SYR  is  also  quicker  than  WRSMOM  for  matrix  fill. 
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This  is  surprising  since  a Galerkin-triangle  moment  method  is  more  compli- 
cated than  pulse  expansion-impulse  weighting.  However,  the  timing  difference 
could  be  attributed  to  differences  in  programming  efficiency. 

Note  that  matrix  fill  is  the  dominating  contributor  to  computer  run 
time  for  N less  than  ^100. 


7.  CONCLUSIONS 

Three  typical  thin-wire  moment  method  computer  programs  ware  compared. 
Two,  OS  and  SYR,  are  basically  Galerkin  with  sinusoids  and  triangles  respec- 
tively. The  third  WRSMOM  employs  pulse  expansion  and  impulse  testing.  These 
programs  were  applied  to  four  thin-wire  antenna  examples,  two  simple  and 
two  complex.  The  parameters  computed  were  driving  point  Impedance  and  far- 
field  radiation  pattern.  The  principal  variable  was  the  number  of  expansion 
functions  N.  The  program  generally  requiring  least  N to  achieve  fair  ('''SdB) 
agreement  with  converged  parameter  values  was  considered  superior.  Thus  the 
main  purpose  was  to  determine  if  a "best"  computer  program  exists  given 
"coarse"  accuracy  as  acceptable  in  predicting  pattern  and  impedance  for 
complex  wire  models. 

The  two  simple  models  (straight,  resonant  wire  and  three-wire  junction) 
were  studied  to  gain  confidence  in  the  programs  and  evaluation  method.  For 
these  examples  OS  was  decidedly  superior  in  computing  driving-point  impedances 
and  all  programs  exhibited  adequate  pattern  convergence  with  SYR  lagging 
a bit. 

The  two  complex  wire  models  were  (a)  an  eighteen  wire  "whip-excited 
twin  fan"  (Figure  19),  and  (b)  a nineteen-wire  "self-excited  twin  fan" 
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(Figure  26).  The  latter  is  a relatively  simple  model  of  a very  complex 
shipboard  antenna.  The  model  includes  an  image  to  simulate  a perfectly  con- 
ducting, flat  ground. 


The  whip-excited  twin  .fan  model  is  representative  of  large  complex 
structures.  It  includes  junctions  formed  by  as  many  as  six  intersecting 
wires.  However,  the  excited  wire  is  relatively  uncomplicated.  For  this 
example  all  three  programs  demonstrated  satisfactory  impedance  convergence 
with  OS  having  a slight  edge.  The  OS  and  WRSMOM  pattern  convergences 
were  about  equal  and  superior  to  that  of  SYR. 

The  self-excited  twin  fan  is  similar  to  the  whip-excited  model  except 
that  the  excited  wires  are  now  on  the  twin  fan  part  of  the  model.  These 
wires  are  quite  short,  'v  0.0025A,  and  they  separate  two  multiple-wire 
Junctions.  First  the  programs  were  applied  directly  to  this  severe 
example.  Only  SYR  provided  reasonable  results.  The  OS  runs  were  repeated 
with  an  optional  time-consuming  matrix  element  computation  in  terms  of 
exponential  Integrals.  This  resulted  in  convergence  of  impedance  and 
pattern  similar  to  that  for  SYR.  The  WRSMOM  results  were  exceptionally 
poor.  However,  this  difficulty  can  often  be  eliminated  by  altering  the 
model.  An  example  of  one  such  "equivalent"  model  was  discussed  in  Section 
6.2. 

The  often  cited  drawback  of  WRSMOM  of  Inadequate  junction  treatment  was 
not  evident  in  these  examples.  A more  plausible  culprit  is  large  adjacent 
segment  length  variations.  This  is  suggested  in  Table  IV  and  Figures  31-33 
where  lengthening  the  feed  wire  on  the  self-excited  twin  fan  model  corrected 
the  WRSMOM  results.  However,  even  adjacent  segment  length  variations  of 


^2:1  brought  reasonable  results. 

One  drawback  of  OS  Is  an  Inability  to  analyze  models  with  wires  of  dif- 
fering radii.  Thus  each  example  exhibited  a uniform  wire  radius.  It  was 
noted,  however,  that  this  constraint  may  prove  minor  in  the  process  of  re- 
ducing a highly  complicated  structure  to  a viable  wire  model. 

The  [Z]  computation  time  for  SYR  was  about  half  that  for  OS  and  even 
WRSMOM  for  the  same  N.  Also  for  the  same  N the  simultaneous  equations  (7) 
solution  time  for  OS  was  about  one-sixch  that  for  SYR  or  WRSMOM.  This  is 
because  SYR  and  WRSMOM  employ  the  slow  matrix  inversion  method.  However, 
quicker  subprograms  to  solve  (7)  can  easily  be  incorporated  in  SYR  and 
WRSMOM. 

The  overall  conclusion  is  that  for  "coarse"  accuracy  in  solving  complex 
wire  models  by  moment  methods,  most  currently  available  programs  are  prob- 
ably comparable.  It  is  important,  however,  that  a convergence  test  be 
performed  for  each  example  regardless  of  the  program. 
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